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. Abstract 

Starting from the lattice realization of the Ising model defined on a strip with integrable 
O I boundary conditions, the exact spectrum (including excited states) of all the local integrals of 

' motion is derived in the continuum limit by means of TBA techniques. It is also possible to 

follow the massive flow of this spectrum between the UV c = 1/2 conformal fixed point and the 
ff^ ■ massive IR theory. The UV expression of the eigenstates of such integrals of motion in terms 

of Virasoro modes is found to have only rational coefficients and their fermionic representation 
turns out to be simply related to the quantum numbers describing the spectrum. 

57' PACS: 11.25.Hf 

^ ! Introduction 

, It is well known that a deep connection exists between integrable models of statistical mechanics 
and integrable quantum field theories. In particular in quantum field theory the Yang Baxter 

■ equation (YBE) plays an important role as a constraint on the 2-particles S*— matrix. On the other 
hand in statistical mechanics the same equation appears as an equation satisfied by the Boltzmann 
Weights. Boltzmann weights satisfying the YBE are then used to build families of commuting 
transfer matrices, which is another quite general feature of integrability for lattice models. 
If we take as a prototype the RSOS models one finds that the Boltzmann Weights can be 
understood as elliptic solutions of the YBE, and actually it is possible to recognize that the solutions 

^ ' one finds are connected with S— matrices in the Sine-Gordon theory [T]. 

■ Furthermore it is well known that upon a suitable restriction of the couplings the Sine-Gordon 
model is equivalent to minimal conformal field theories [2] Oil]) and actually also the A„ RSOS 
have been shown to be in the universality class of minimal CFTs [5]. 

Nonetheless the A„ are in correspondence with CFTs only in an appropriate continuum limit 
called UV scaling limit. In general there will be a continuum scaling limit depending on a mass 
parameter /i which will generate a RG flow to a massive IR theory, where the relevant processes 
will essentially be the scattering of kinks. 

The continuum field theory corresponding to A„, models can be interpreted as a (f)i,3 thermal 
perturbation of the Mn.n+i minimal conformal field theory [3] . Such a perturbation is known to 
be integrable [6] [7], this means that there exists an infinite number of commuting currents which 
remain conserved in the perturbed theory. 
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In particular the first conserved quantity is the energy, if one considers its value on the vacuum 
state it is well known that this is proportional to the central charge of the underlying CFT in the 
UV limit. The flow to the IR of such a quantity represents an example of the famous c— theorem. 
A powerful tool for having access to the vacuum energy is the Thermodynamic Bethe Ansatz 
(TBA)(see for example [8]), which in some cases has been generalized to excited states [5]|9] [10j|llj . 
In this work we derive the excited TEA equations for the A3 model on a strip with integrable 
boundary conditions by diagonalizing the transfer matrix. We then proceed to define the continuum 
scaling limit of the transfer matrix eigenvalues which we then use as generating functions for some 
quantities which we eventually identify with the conserved quantities of the thermally perturbed 
conformal field theory. 

It is then possible to analytically follow all the conserved quantities along the massive flow to the 
IR theory. Comparison of the results which are obtained in the UV limit with the spectrum of 
the BLZ local integrals of motion provides an exact identification of the conserved quantities and 
allows to put the lattice boundary conditions in correspondence with the CFT operator content of 
the theory. 

The eigenstates of the BLZ integrals of motion are computed in their Virasoro form, and once 
expressed in terms of the fermion field turn out to be labelled by the same quantum numbers 
which label the exact formula for their eigenvalues which has been derived independently through 
TBA. 



1 The ^3 Model 

The A3 model is a lattice model which provides a convenient realization of the Ising model. It 
is built on a square lattice where to each site j is assigned a height variable aj £ {1,2,3}. The 
local height variables aj are constrained to satisfy an adjacency rule which holds for i,j nearest 
neighbours: 

\ai — aj\ = 1 (1-1) 



The model is characterized by its Boltzmann Weights [12], which can be of two different types: 
bulk weights and boundary weights. The only nonvanishing bulk weights are: 
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where A = 7r/4 is the so cahed crossing parameter, u is the spectral parameter and the ^ are related 

to the choice of boundary condition. 

We also define the elliptic theta functions of nome q as: 



h{u,q) 
h{u,q) 
h{u,q) 



V^^ sin((2A: + l)n) \q\ < 1 

fc=0 

oo 

2gi/4^gfc(fc+i)cos((2A: + l) u) \q\ < 1 

k=0 

oo 

l + 2^(-l)'=g'='cos(2A:n) 



(1.5) 



< 1 



k=l 



this is the so called g-series, which will prove more useful to our goal, more typical definitions of 
these functions are given in terms of infinite products. 

It is important to remark that the role of the nome q is to control the criticality of the model, 
which becomes critical as g — > 0. In what will follow we will focus on the region < q < 1 which 
is the so called regime III of [5], actually this regime corresponds to a low temperature phase, but 
because of the duality symmetry of the Ising model this is the same as a high temperature phase. 
Now, in terms of the above objects it is known that the model admits a transfer matrix description 
(on a lattice of width N): 
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(1.6) 
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such a transfer matrices form a one parameter commuting family with respect to the spectral 
parameter u, and it is well known that this property makes the model integrable. 
The transfer matrix T(n) satisfies the following functional equation: 



T{u)T{u + A) = 1 + diu) FNiu) Siu, Cl,^r) = 



(1.7) 
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, , , 6ii(2'u-2A)6'i(2ii + 2A) , , ^ , , , , , 

^{u,U,^R) = -^-77, TTT-TT^ — 9, 4l, «L)-^fl(^^, 9, (1-9) 



6'i(2'u- A)6'i(2'u + A) 
where d is a matrix proportional to the identity that takes the form: 

27V 

d(u,g) = l(-l)^' 



v{u)Qi(u - 2A) 



y(u-X)dx{u^X) 



L(2n) 
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(1.10) 



being 



OAju + 7r/4 - ^iWiju - 7r/4 - Cl)04{u + 7r/4 + Cl)Oa{u - 7r/4 + ^l) 

04{u - aL7r/4 - ^l)04(.u + ultt/A + a)^4(n - (ol + 2)7r/4 - a)^4(u + (ol - 2)7r/4 + Cl) 

(1.11) 
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^ + 7r/4 - ^R)ei{u - 7r/4 - CnW^ju + 7r/4 + Cr)94{u - tt/A + Cr) 

^ 6'4(u - aRTT/A - iR)ei{u + a^j7r/4 + iR)ei{u - {aR + 2)7r/4 - ^R)e4{u + {ur - 2)tt/A + ^r) 

(1.12) 

AL[u,q,^L,aL) = e4{X)^ ^ 

- / ^ s 04{u-CR)e4{u + ^R)e4{u + aRTT/A + CRW4{u-aR7r/A-^R) 
■Ar{u, q, (r, aR,) = e^(Xj^ ^ 

The phases in the B terms may seem strange, but they turn out to be necessary. This fact has 

been observed also in [13] where the TBA equations for the A3 model were derived. 

Such a matrix satisfies a functional equation which for obvious reasons is called the Inversion 

Equation: 

d{u)d{u + X) = l (1.15) 

As a consequence of the simple form of the d matrix (which for more complicated models is not 
diagonal but is expressed in terms of T itself), we have that the Functional equation written in 
terms of the eigenvalues T of T is independent of the eigenvalue under consideration. 
Before moving on to discuss the TBA equations, it is useful to spend some word to comment on 
the periodicities of the transfer matrix T{u). Such periodicities come directly from the properties 
of the elliptic 9 functions and read: 

T(n + 7r) = T(ti) (1.16) 

T{u-ilogq) =T{u) (1.17) 

As a consequence we have that T is a doubly periodic function which is completely defined by its 
analytic properties inside a rectangle that we may take as: 

(-^,^A)xi(ilog(g),-ilog(g)) (1.18) 

If we now consider the functional equation ()1.7p it is clear that Z?'s periodicities are inherited by the 
righthand side J-'{u) so that the object of our interest will be the zeroes of J-' inside the periodicity 
rectangle (ll.lSj) . 

Such zeroes can be shown by numerical analisys (and indeed analytically in the critical limit [H] 
) to be organized on lines parallel to n = A/2 + ix with periodicity A. It can be argued (see again 
[H] for the critical case) that as a consequence of the periodicities and of the structure of the 
functional equation (|1.7|) the eigenvalues T(n) must have zeroes which are organized in structures 
called 1 — strings and 2 — strings. 

1-strings are just single zeroes of real part A/2 and imaginary part < < —1/2 log g such that: 

T{\/2±ivk) = k = l,...,m (1.19) 

where m denotes the number of 1-strings, while 2 strings are couples of zeroes sharing the same 
imaginary part, while their real part takes the values A/2 it A, and we shall call their number n. 

It is clear that a 1-string is a zero both for T{u) and for J-{u), therefore is will be convenient 
in our case to look for such zeroes in the expression for T. Such zeroes coul in principle be found in 
either one of the three factors of which T is composed, but since (as we shall shortly see) i^jv and 
S are going to be eliminated in the scaling limit 1-strings will essentially be zeroes of the (1 -|- d{u)) 
term. It is finally worth mentioning that for finite it is possible to give a characterization of 
the states (i.e. transfer matrix eigenvalues) in terms of a non increasing sequence of numbers 
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Figure 1. Example of the structure of zeroes labeled by the topological number {0,0} 



{Ii, I2, ■ ■ ■ , Im} called quantum numbers which express the position of 1-strings related to the 
position of 2-strings. 

Each /fc tells us how many 2-strings the k-th 1-string has to go through in order to reach its position 
in the pattern of zeroes starting from the configuration where all the 2-strings are heaped on the 
bottom, so that ordering the imaginary part of the 1-strings v/^ into an increasing sequence {vk}^^i 
we have that the quantum numbers {Ik} must necessarily arrange into a non increasing succession. 
Clearly, the 1^ have to satisfy the following constraint: 



Such a characterization of the eigenvalues in terms of 1-strings and 2-strings also happens to give 
us a natural criterion for ordering the states, first of all we order the states by their increasing 
m value, the ordering between equal m states is done so that the state with all the 2-strings at 
the bottom of the tower comes first, and then each time a 2-string is pushed over a 1-string the 
"energy" increases by one "unit" . 

We shall see that in the continuum limit a more natural set of quantum numbers will arise to 
describe the pattern of zeroes. 



In this section we are going to derive the excited state TBA equations for the model by solving 
(jl.7p . Considerable work has been done in the past on the excited states TBA, here we will 
essentially follow the work of [5]. 

First of all let us recall the form of the Functional equation (jl.7p . we then define an x coordinate 
in the following way: 



(1.20) 



2 Excited TBA Equations 





(2.1) 



where we are going to solve (jl.7p for the following values of x: 



X G (21ogg, — 21ogg) 



(2.2) 



for convenience we will rewrite (jl.7|) after applying a traslation: 



u^u-\: T{u-^)T{u + ^)=Hu-\) (2.3) 
we then use (I2.ip to write ()1.7p in the following form: 

''=^ + '1' T,{x + i^-)T,{x-t^-)=nn-\)=:F{i-/^Tr{x) (2.4) 

At this point we could be tempted to follow the solution method used in [5] and try to Fourier- 
expand the logaritmic derivative of (|2.4|) . anyway before being allowed to do so, we have to remove 
the zeroes of Ti(x) in order to deal with an analytic function for which a Fourier expansion does 
make sense. 

Now, if we consider what has been said in the previous section about the position of the zeroes, 
one observes that for real x, |x| < — 21ogg, the function Ti[x) which are due to the presence of 
1-strings. 

In order to reach our result we define the function 

(2.5) 

We observe that the p function satisfies an equation which is similar to Ti(x): 

vr vr 

p[x + i-,Vk) p{x-i-,Vk) = l (2.6) 

furthermore we observe that p can be used to collect all the zeroes (for real x) of Ti through the 
product: 

m 

^p{x,Vk)p{x,-Vk) (2.7) 
fc=i 

so that we can assert that the function 

T.M^) = J—— — - = n -4..),,=) ii(i(rTw) 

YYp{x,Vk)p{x,-Vk) 

k=l 

does not have zeroes for real x \x\ < — 21ogg (ANZ stands for analytic and not zero). 

We then observe that as a consequence of ()2.6p T4 )\fz still satisfies the Functional equation: 

rANz(x + i|) Tanz(x - i|) = J'lix) (2.9) 

so that now one is authorized to fourier-expand the logarithmic derivative of the above equa- 
tion. Bymeans of some algebra one can determine Tanz{x), and thus Ti(x) to be: 

m 

logTi{x) = Y,^og[pix,Vk)p{x,-Vk)]+k*logJ^i + D (2.10) 

k=l 

Where k{x) is a convolution kernel defined as: 



1 — A e 2 log? 

k(x-y) = -— V 5 ^ (2.11) 

41oe:o ^ __fe2£_ fc^^a v j 

^^k = -O0e 41ogg _)_ g41og5 
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The convolution kernel k{x) can be computed in terms of Elliptic 6 functions. It has been computed 
in [5] to have the following form: 



2tt02{ix, q^) 



Finally if we recall J-{u)^s definition we can write: 

\ogJ^,{x) = logJ^ii^) = log (l + + logF^iij) + log5(i|,a,eR) (2.13) 

so that remembering the expression for Fj\f{ij) we observe that its logarithm gives a contribution 
proportional to N and thus can be identified with a Bulk Energy term (which diverges in the 
thermodynamic limit N oo), while S{ij,£,L,^R) gives a contribution independent of N which 
anyway diverges as we approach the critical regime {q — > 0) and can thus be identified with a 
Surface Energy. 

Now, since in the next section we are going to deal with the continuum limit of log T(ii) (which 
consists both of — > oo and q ^ 0), it is natural to conclude by defining a subtracted Energy so 
as to give rise to meaningful quantities in the continuum limit: 

logTfi,ite(x) = logTi(x) - A; * log [F;v(i|)5(i|, a, e^)] (2.14) 

and, more explicitly: 

m 

logrfinite(a;) = ^log[p(x,ffc)p(x, -Vk)] + A; * log (l + d{i-)) (2.15) 



Where we observe that the constant of integration D has finally disappeared in the subtraction, 
corresponding to a shift in the vacuum energy. 

3 Scaling Limit 

Taking a scaling limit of a lattice model essentially means considering its critical behaviour in the 
thermodynamical limit. 

such a double limit {N — > oo, ^ — > 0) can in principle be computed along infinite paths, however 
it can be shown [5] that there exists a prescription which allows us to obtain a meaningful limit 

for log D finite ■ 

such a prescription has the form: 

q = t^ (3.1) 

u = ^{x + logN) (3.2) 
where the size and the reduced temperature t satisfy: 

iV ^ oo (3.3) 

Nt" = fi (3.4) 

1/ is understood as the critical exponent of the correlation length, which for the A„ models in 
regime III of [5] is known to be: 

= —— (3.5) 
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The II parameter plays the role a regulator for the continuum system, and it will be used to generate 

a massive RG flow connecting the UV {fi = 0) and IR = +00) fixed points. 

Such a regulator can be thought of as arising from the product of a mass m and a length R: 



being 



IJ = -mR (3-6) 
R= lim Nl (3.7) 

t" 

m = lim — (3.8) 
and / is understood as the lattice spacing. 

Our first goal is to build a continuum version of d(n, q) which we will use to expand the continuum 
transfer matrix as a series whose coefficients will be the integrals of motion, we will then use the 
information we will gain to dicuss the full continuum expression for \ogT^aitc{x) ■ 
The limit we are essentially interested in computing is the following: 

d{x) :=^lim^d(l(x + logiV),(^)^) (3.9) 

From now on we shall fix the right boundary to a/j = 1 so that already before going into the scaling 
limit the R boundary term is fixed to: 

Br = -1 (3.10) 

Therefore from now on the only surviving boundary height will be simply called a and we shall 
realize in what will follow that a = s where s is the Kac label of A^^s, and of course in this case 
r = aR = l. 

In order to achieve the correct scaling behaviour for the boundary term one has to postulate the 
following scaling behaviour for the boundary parameter ^ (we are dropping the subscript L): 

e~r + ^iogiv (3.11) 

It is also useful to split the boundary term as: 

B = e^™e+S- (3.12) 

so that B'^ depends only on u + ^ and similarly B~ depends only on n — ^. Proceeding further one 
finds: 

^+ _ 1 - /i^ete~^'^^*+^) 1 - ;u^ete~^^(^*~^^ 

" 1 _ ^^ef e-2^«*+'^f ) 1 - ^lef e-2^«*-(«+2)f ) ^ ' ^ 

- 1 - uiete~^*(^~^*) 1 - /iiete^*(?+^*^ ,^ , 

B~ = — — (3 14) 

l_/,lefe2^(«f+e) i_^lefe2^((«+2)|+e) ' 

One then notices that in the UV limit /U ^ the boundary B flows into 

B ~ e^™ (3.15) 

and takes the value —1 for a = 1, 3 and +1 for a = 2. 

After carrying out the calculation one finds the following continuum version of d : 

d{x, /i) = e-S/^-^M^+iog/.) tanh^ (^^±Mi!) {^BlBr) (3.16) 
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One should at this point discuss the scahng Umit of the convolution term in log Tgnitcl^^)) such a 
calculation turns out to yield: 

lim k*log(l + dr-{x + logN))) = r ^ -\og{l +d{y))dy (3.17) 

N^oo,q^O \ \4 // ZTT COSIl[X — y) 

We can now deal with the scaling limit of the excitations. 

As we approach the scaling limit the 1-strings will have the following asymptotic behaviour (yk is 
the finite part): 

Avk^yk + logN (3.18) 

The 1-string term will then become: 

±iogipix,v,)pix, -V,)) ~ flog (tanh (i^^I^i±i^^) tanh ( " + ^^ + ^°g^ )) (3.19) 

k=l k=l 

4 Expansion 

4.1 Study of the Ground State 

We first of all begin by studying the behaviour of the continuum ground state eigenvalue of the 
transfer matrix; in the previous section we have shown it takes the following form: 

logDM=pi?^J2i%tM (4.1) 
* ^ ' y_„ 2iT C08h(i -y) ^ ' 



such an expression, following the spirit of [15] has to be expanded in the following series: 

oo 

log D{x) = -J2 Cnhn-i (/i)e(2"-i)- (4.2) 

n=l 

which in our case yields the following expression for the vacuum integrals of motion: 

(_-\\n r+oo 

CnlZ-i{^^) = / dye-<^'^-''>y log(l + d{y)) (4.3) 

we're now going to study this expression in the UV and IR limits. 

We start by observing that the above expression can be manipulated into the form: 



^ ,.ac _ (-i)v^ V r ^S(2"-i)p-(*+^)r^^fVR B 



k 



(4.4) 



Where we once we fix the right boundary to = 1 Br disappears from the equations. 

One at this point decides to get rid of ^* fixing it to zero. Actually one could keep it, and let it 

scale once again as /i ^ so that, for a suitable A: 

^Ag4ir = ^ (4.5) 

The parameter r would then generate a flow between different conformal boundary conditions, see 
for example |16] . 
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Now, by means of standard techniques one can prove that the following inequality holds for all 
values of a: 

2i,(2"-i) ^ 1 

\CnlZ-i{f^)\ < -^^-^ E l^l-2n{^k^^) (4.6) 

k=l 

where the Ki{z) are the modified Bessel functions of the second kind. A study of the large /i 
asymptotics of the above series allows one to conclude that the ground states of the integrals of 
motion decay exponentially in the IR limit. 

We are now ready to move our attention to the UV asymptotic behaviour. 
In the limit /i — > it is not difficult to show that: 

Cn/2T-i(A^) ~ ^^^r(2n - l)Li2„(e-'^) (4.7) 

where 

oo ^ 
k=l 

is the Polylogarithm function. It is worth spending a word to observe that considering the energy 
Ii one finds a dilogarithm of a phase, actually it is well known that the central charge is usually 
proportional to a dilogarithm. Furthermore one notices that all the vacuum expectation values 
of the integrals of motion are proportional to polylogarithms, this seems to be a rather general 
structure. One could ask himself if such polylogarithms satisfy sum rules similar to those holding 
for dilogarithms. 



4.2 Excited States 

We now ask ourselves what is the behaviour of the excitation terms in the UV and IR and compare 

it with the ground states. 

Now we want to expand the 1-string term. 

The expansion is readily obtained if one considers the following identity: 

[ -\ f\ °° f2k-l 



and for example writes: 



logtanh ^ ^fc + ^°g/^ =^^_2^e(^"-^)"- ^ (4.10) 

n=l 

Clearly one obtains the following result: 

1 - String) = - E _ ( E (^"^'""'^^' + ^2(2n-i)e(2n-i)s/.)^ ^4,11) 

n=l k=l 

In order to proceed further it is necessary to have a closer look at the asymptotic behaviour of the 
Vk- 

Let us recall the equation satisfied by the yk'- 

d{yk-il)=-l (4.12) 
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so that by taking the logarithm of both sides, using the expression for d given previously we have: 

- 4fi{e-^'\y^fi - e''^e-y^^r') + 21og ^ ^2iX^.i^.yJ + log(^L^fl) = ivrn^ (4.13) 

where the are odd numbers. 

We now introduce the following function: 

gkifJ.) := fiey^ (4.14) 

we are now interested in expressing the yk equation in terms of this new function g^, in order to 
do so we first rewrite the boundary term as: 

Hvk - 4) = %'.^!]\ , (4.15) 

2 [gk + - ^i9k[l + cos(7ra)) 

so that the yk equation becomes simpler and reads: 

4i^M:il) + log ( — r^) = + 1 - a) (4.16) 

9k \{gk + ty -2tgk{l + cos{TTa))/ 

this equation gives us the inverse function fj-{gk), and we are interested in its behaviour as ^ ^ Oand 
/i +00. In order to reach our goal it suffices to pick the branch of the fi{gk) function which 
passes through the origin, gfc — > corresponds to the UV limit, whereas 5^ — > 1 is the IR limit. 
Expanding fi around g^ = we obtain: 

fir^jia-l-nk)gk + 0{gl) (4.17) 

If we now decide to rewrite the 1-string term as follows: 

00 „ I2n~l)x -I 

(1 - strin,) = - Y: ^T^M-"-" E ii"-" + T^Tt) ("8) 

n=l k=l 9k 

we conclude by applying the above UV expansion that the expression 

9 k ' 

is UV limited, so that in this limit the excitations have the same scaling behaviour as the ground 
state. 

On the other hand the IR excitations cannot avoid to grow faster than the ground state term, this 
observation united to the fact that gk ^ 1 dictates the particular structure of the IR spectrum. 



4.3 Full Expansion 

We finally have arrived at the point of writing the analytic expression for all the integrals of motion 
of the model, such expression reads: 

2^(2.-1) ^ , 1 



c„/,„-.(.)=^E(r-" + T^)+ 



fc=i 9k 

+ 



(4.20) 

/ -i\n r+00 



vr 
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In the limit fi ^ oo, — > 1 the ground state drops off exponentiahy, so that considering only the 
excitations it is immediate to realize that: 

C„/2n-l(^) ~ -^^—r (4.21) 

this spectrum happens to be integrally spaced, and looses all memory of the quantum numbers 
aside from the length of the sequence of the Uk- 

It is worth remarking that this result should not surprise us very much since it is very similar 
to what has been achieved in [T7| for the A4 model corresponding to the tricritical Ising Model 
universality class. 

In the UV limit it is clear from what we said so far that the integrals of motion have the following 
behaviour: 

fc=i (4.22) 
+ (_l)"e^™((l - e^™)2-2" - l)r(2n - 1) ^^^"^ 



^2n 



For a = 1 the energy takes the following form 

1 1 

Ci/i(A^)~vr(-J](-nfc)--) (4.23) 
k=i 



For a = 2 we have: 



For a = 3 we have: 



Cihif^) ~ ^(^ f;(l - Uk) + ^) (4.24) 

k=l 



1 1 

Ci/i(M)~vr(-E(2-nfc)--) (4.25) 



2 48 
fc=i 

It is readily recognized that these formulae are in agreement respectively with the h = 0,h = 1/16 
and h = 1/2 sectors of the minimal model A^3,4 if we choose: 

Ci = vr (4.26) 

One notices that in the vacuum sector m must be even whereas in the 1/2 sector m must be odd 

for trivial reasons. In the 1/16 sector m must be odd but this fact is less trivial to understand 

from the formula, let us simply say that m— parity is fixed in the sector and is odd because for the 

highest weight state the only possible quantum number must satisfy 1 — ni = 0. 

We will understand better the structure of the quantum numbers in the next section. 

If we now want to compute the constants C2 , C3 in the vacuum sector of the model this can be 

done by using the explicit expressions for the integrals of motion, which can be found in |18j : 



Ii = ^0 - ^ (4.27) 
I3 = 2 f; L„„L„ + Ll- ^Lo + (4.28) 

n=l 
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3 °° 

I5 = Sm+n+p, : LmLnLp : +— Li_2ri-f'2n-l + 

m,n,pS:Z n=l 



+ V - ^ - 1 ^nL_„ - + (c + 2)(3c + 20) (4.29) 

^\^6 4 y"" 8° 576 

_ c(3c + 14)(7c + 68) 
290304 

where the : : denotes Conformal Normal Ordering which can be obtained by arranging all the 

in an increasing sequence with respect to n. 

So that we have the following vacuum expectation values: 

49 

(0130) = — — (4.30) 
\ Ml / ;l1520 ^ ' 

/0|l5|0) = (4.31) 

\ \ ^\ ' 2322432 ^ ' 

which gives us the following equations: 

•^^CI'^IO).-^™^ ,4,33) 



SO that we readily get: 



C2 = - (4.34) 
Cs = (4.35) 



This values for the vacuum constants C„ actually can be extracted from [19], so that in general we 
have for the vacuum sector the following expression: 

a=^"^^""--^""^(^--^) (4.36) 
n!r(3n-i) ^ ^ 

Actually these values of the C„ are computed from the vacuum sector, but direct calculation allows 
one to verify that they are independent of the sector. 

We stress that the above formulas describe exactly the conformal spectrum of the minimal model 



5 Fermionic modes and TBA quantum numbers 

Actually the quantum numbers themselves have a very simple interpretation in terms of 
fermionic modes. 

To understand this one needs only to remember that the stress energy tensor for the Ising model 
is built out of the fermion field as: 

T{z) := i : ^{z)di^{z) : (5.1) 
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CFT State Fermionic State TBA State 


0) 0) 





2L_2 0) i/'_3i/'_i 


0) (-1,-3) 




0) (-1,-5) 


1 1 T T T 1 f * *-| 

fL_4 0) + fL?.2|0) ^-l^-l 


0) (-3,-5) 
0) (-1,-7) 


|i-5|0) - |L-3i-2|0) l/'_lV'-3 

fL_5|0) + fL-3L-2|0) V'-aV'-i 


0) (-3,-7) 
0) (-1,-9) 


f^L_6|0) + fL-4i-2|0) - §LU\0} V'-iV'- j 
iL-6|0) - iL_4L-.2|0) + j^L-ajO) V'-fV'-f 
^L_6|0) + ^L-4L-2|0) + ^Ll,\0) V-4V'4 


0) (-5,-7) 
0) (-3,-9) 
|0) (-1,-11) 



Table 1. U.V. state correspondence CFT — >TBA for the h = sector 





CFT State 






Fermionic State 


TBA State 




1/2) 


0) 


(1) 




L_ 


-i|l/2) 










3 |0) 


(-1) 






-2|l/2) 








i,_ 


5 |0) 


(-3) 




kL 


-3|l/2) 








^_ 


7 0) 
2 1 ' 


(-5) 




.4jl/2) + iL_3i- 
-4|l/2)-|L_3i- 


i|l/2) 
i|l/2) 




i'- 


V'-||o) 

5ih 3ih i\0) 


(-7) 
(1,-1,-3) 




-5|l/2) 
-5|l/2) 


+ ^L-aL. 

— |l/_4L_ 


i|l/2) 
i!l/2) 






2 


11 |o) 

lip ilO) 

2 2 1 ' 


(-9) 
(1,-1,-5) 


iL-6|l/2) + 
\L-,\1I2)- 
|i-6|l/2) + 




-i|l/2) + 
-i|l/2) + 
-i|l/2)- 


^L-aL- 


-2|l/2) 
-2|l/2) 
-2|l/2) 




!/> 13 |0) 

9 1/) 3tb 1 0) 
2 2 2 ' 

7 5 1 0) 
2 2 2 1 ' 


(-11) 
(1,-1,-7) 
(1,-3,-5) 



Table 2. U.V. state correspondence CFT — >TBA for the h = 1/2 sector 



CFT State 



Fermionic State 



TBA State 



1/16>_ 



2V2L-i|l/16) 



V2i>o 1/16)_ 



(1) 



^L-2|l/16)_ 



^-l|l/16). 



(-1) 



fL_3|l/16)_ - i^L-2L-ijl/16)__ 
2^L_3|1/16)_ + ^L_2L-i|l/16)_ 



^-2|1/16)^ (-3) 
V2V'-2V'-iV'o|l/16)^ (1,-1,-3) 
V-3|l/16>^ (-5) 



V2t/>_3^A-i^o|l/16)^ (1,-1,-5) 
^-4 1 1/16), (-7) 



fL-4|l/16)_ -L-3L-i|l/16)_ 
^L-4|l/16)_ + -^L-3L-i\l/16)_ 



8%/2 



2 

28-^ 



iL_5|l/16)_ + |L_4L-i|l/16)_ - L-3L-2il/16)_ 
|1/16)_ - §L-aL.i\1/16)_ + j|L_3L_2|l/16)_ 
L-4L-i|l/16)_ + ^L-3i-2|l/16) 



^/.-5|1/16)_ -f^i.-4i.-l|i/iD;_ +^^-3^-2|l/iD; 

^L_6|l/16)_ - ^L.5L-i|l/16)_ + -i5L_4L_;|l/16)_ - ^^Li\l/16)_ 
i|L.6|l/16)_ + ^L_5i-i|l/16>_ - AL_,L_,|i/i6)_ _ ^Li3|i/16)_ 

;|L_6|l/16)_ - |L_5L-l|l/16)_ + |L?.3|l/16)_ 

44f^^-^|Vl6)- + sBf ^-^^-^IVl6)- + 7^^-4^-2|l/16)_ + Titr2L'-s\l/lG)_ 



^/2t/>_3^/;-2^o|l/16)^ (1,-3,-5) 
V2t/'-4^/'-i^o|l/16)^ (1,-1,-7) 
^-5|1/16), (-9) 



^_3^-2V-i|l/16)+ (-1,-3,-5) 

V2V'-4V'-2^o|l/16)^ (1,-3,-7) 

^/2V-5^/'-l^o|l/16)^ (1,-1,-9) 
V-6|l/16), (-11) 



Table 3. U.V. state correspondence CFT — >TBA for the h = 1/16 sector 
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TBA State Eigenvalue Is Eigenvalue 
59 -4433 







-1,-3) 



-1,-5) 



-3,-5) 



-1,-7) 



-1,-9) 



-3,-7) 



-5,-7) 



11520 2322432 

4Y089 17581135 

11520 2322432 

211729 225292639 

11520 2322432 

2554U9 242734063 

11520 2322432 

577969 1211381743 

11520 2322432 

1226449 4255847167 

11520 2322432 

621649 1228823167 

11520 2322432 

786289 1436534671 



-3,-9) 



(-1,-11) 



- rMfe 

11520 2322432 

2237809 11607335311" 

11520 2322432 



Table 4. Table of eigenvalues in the vacuum sector 



TBA State I3 Eigenvalue Eigenvalue 
TTW 67639 



(1) 



(-1) 



(-3) 



-5) 



11520 2322432 

45409 17509063 

11520 2322432 

21()049 225220567 

11520 2322432 



(-7) 



576289 1211309671 



(1,-1,-3) 



(-9) 



(1,-1,-5) 



(-11) 



(1,-1,-7) 



(1,-3,-5) 



11520 2322432 

257089 242806135 

11520 2322432 

2236129 11607263239" 

11520 2322432 

623329 1228895239 " 

11520 2 322432 

3691009 26759824663 ' 

11,520 2322432 

1271809 4273361)663 

11520 2322432 

787969 143660674 



787969 1436606743 

11520 2322432 



Table 5. Table of eigenvalues in the 1/2 sector 



TBA State 


I3 Eigenvalue 


Is Eigenvalue 


(1) 


-7 
1440 


143 

— m — 


(-1) 


1673 
1440 


72576 


(-3) 


13433 
1440 


2306447 
72576 


(1,-1,-3) 


15113 
1440 


2378519 
72576 


(-5) 


45353 
1440 


17*513639 
72576 


(1,-1,-5) 


47033 
1440 


17*585711 
72576 


(-7) 


107513 


73801871 


(1,-3,-5) 


— m — 

1440 


nMfe 

72576 


(1,-1,-7) 


109193 
1440 


73873943 


(-9) 


209993 
1440 


72576 


(-1,-3,-5) 


60473 
1440 


19892015 
72576 


(1,-3,-7) 


120953 
1440 


76108175 
72576 


(1,-1,-9) 


211673 
1440 


225297215 
72576 


(-11) 


362873 
1440 


560432015 
72576 



Table 6. Table of eigenvalues in the 1/16 sector 



15 



So that by introducing the well known mode expansion 

= (5.2) 

one gets 

-C-n = + -) : V'n-fcV'fc : (5.3) 

k 

Where the fermionic modes will have a half integer index when we will be working in the 0, 1/2 
sectors, whereas the index will be integer in the twisted 1/16 sector. 

It is then just a matter of unraveling the normal ordering and using the fermionic algebra {ipn, V'm} = 
Sn+m,o to work out the fermionic expression for the eigenvectors of the integrals of motion. 
If we consider for example the sixth level of descendance in the vacuum sector we have: 

20L„6|0) + 24L„4L_2|0) + 5L23I0) = 112'i/'_ n V'_ 1 10) (5.4) 

4L_6|0) - 8L_4L_2|0) + ^^.sjO) = 16V'_|V'_| |0> (5.5) 

20L„6|0) + 24L_4L_2|0) - 23L2,3|0> = 56?/'_7 V^^i |0) (5.6) 

So that comparing with tabled] we see that ^ are simply the labels of the fermionic modes, and we 
can easily either guess or explicitly work out straightforwardly the form of all the other eigenstates 
which is obvious aside from a normalization. 
Similarly in the 1/2 sector one has at level 6: 

36L_6|l/2> + 27L_5L_i|l/2) +8L_4L_2|l/2> = 288Vj„ |0> (5.7) 

36L_6|l/2> -45L_5L_i|l/2) + 8L_4-^^-2|l/2> = 144^_9 -^.s V'_ 1 |o) (5.8) 

108L_.6|l/2) + 81L_5L_i|l/2> - 26L_4L_2|l/2) = 72?/'_7 V'™ 1 |o) (5.9) 

so that the fermionic modes have indexes which are simply "''^^ and one understands how the 
fermionic represetation of table [2] should be. 

The 1/16 sector has to be worked out expanding the Virasoro modes in integer fermionic modes, 
this is not very different from the previous situation, except for the presence of the zero-mode ipQ 
which generates the zero mode algebra: 

^o|l/16>_ = -^|l/16>^ (5.10) 

V'o|l/16>^ = -^|l/16>_ (5.11) 

So that actually there are 2 1/16 vacua, one has no 1-strings (|l/16)^) and the other has 1 1-string 

(|l/16)_). The fermionic modes have indexes which are simply "''^"^ , so that the (1) quantum 
number is recognized as coming from the insertion of ipQ. 
For example one has: 

|l/16>_ = y2?^o|l/16>^ (5.12) 
L_i|l/16>_ = ^V'-i|l/16>+ (5.13) 
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L_2|1/16>_ = ^V^_2|1/16>^ (5.14) 

So that comparing with table [3] we easily understand the how things work out and, once again, 
aside from the normalizations one knows perfectly from the beginning which results he will find 
upon expressing eigenstates in terms of fermionic modes. 

It has to be remarked that the fermionic description of the eigenstates of the RSOS transfer matrix 
fits naturally into the description of [SOj, where it is shown that the TBA quantum numbers can 
be described also in terms of "fermionic paths" . 

Finally one notices that the TBA quantum numbers appearing in the 1/16 and in the 1/2 sector 
are actually the same, this can be understood from the existence of a boundary flow connecting 
the 2 sectors, see for example [21] . 



6 Conclusion 

In this work we have derived the excited TBA equations for the ^3 model defined on a strip with 
integrable boundary conditions. The continuum limit of the TBA equations was then derived and 
an axpansion defined which led to integrals of motion. Such integrals of motion were identified to 
be the BLZ local integrals of motion. The identification was carried out by comparing the exact 
diagonalization of the Virasoro expressions with the spectrum derived form TBA. The eigenvalues 
were found to be exactly the same and the eigenvectors, once expressed in terms of fermionic 
modes, turned out to be labelled by the same quantum numbers as the eigenvalues. 
It has to be remarked that this is the first time that an exact diagonalization of all the BLZ local 
integrals of motion has been carried out in a particular case including all the excited states. 
A more detailed analisys of the spectra obtained would surely be interesting and will be the subject 
of successive work. 
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